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عاطقباله۷ عصتصتهاط9اه ما وهع1 (800) مممایعه متصهع تن آزمو مه باتصتلعی آنمو ام جمت‌طاتاوتل لهتاهموه مطا 0۶ 1۱0۱۲۱6026 
5 20 ۱۷۵۲۱0۶ عط 0۶ 0عنط) ه مقطا میم .عمتاز نامه لهعبانهانمه 10۲ عومنلوصص ممزوتمع مصز 7۵ تامع1ه وا عقطا مت1010۳۲۵ 
ماطهمتهامده 0۶ اصمصصمماهبعل مط) فاصههج قمع رومم ۵ ممتامتالمتن صرح همع مط فصمهم‌نطا مامتها ملقی رها 2100160 
-تصصمو مه دنه ما ومعصعمالقده کصهمگنصونه ومومم ولژمو صا کصمنجمی ۶0 از تامنلجمی لهمزتامعاه طعنط م1 ,عتااتاهتتهه 
صانلمو ۵ وزم۷ع۱ طونط کزطانن ماه والمو منلمو مه مصتلهو .صمتاعملمم تحااتهاوه مصتامهم‌صا رامعنع روجملعع۲ 2710 
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صح ره 1 ,تفه نصا ۸220 متصعاوا بطعصهرظ طمنوعوعع همه ممعصمتهگ راصهص‌تهمه1۳ ممممتم6 9011 مطا 0۶ 550۲ع]۳۲۵ -3 
00 ۲۸۲۲۵۱۹10۶ مه ممتاه»۱0۲ ردم‌نهعومک امتانه‌تنهه ,ماتاتافص؟ دعوم تمه ۷۷ مه 5011 ,تمووع]۲۳0 ]ممادزدوظ -4 
صح ,رمک ,رمتتزض) 

+6721 رواتوته۲تص ۸220 متصصهاق بطم‌صهرظ موم مه عمصملمو: رکصمصاندمونا معصمتمو 501 عطا ۵۶ 950۲ع]۲۳۵ ]صداوزوع۸ -5 
۱۹ 

10]: 8 


۶ نشریه آب و خاک. جلد ۰۳۷ شماره ۲. خرداد - تیر ۱۴۰۲ 
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۰ 13 ما اممجصم‌همصممط 5۵11 م] ۲۵۱2/60 
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نقشه‌برداری رقومی تغییرات سطحی و زیرسطحی کربن آلی و شوری خاک در بخشی از اراضی 
دشت قزوین (مطالعه موردی: مناطق آبیک و نظر آباد) 


گردآفرین رضایی (- فریدون سرمدیان "*- علی محمدی تر کاشوند "- حواد سید محمدی *- مریم مرعشی علی آبادی ٩‏ 
تاریخ دریافت: ۱۴۰۱/۱۰/۱۲ 
تاریخ بازنگری: ۱۴۰۱/۱۱/۱۷ 


تاریخ پذیرش: ۱۴۰۱/۱۱/۲۴ 


چکیده 


آگاهی از توزیع مکانی شوری و کربن آلی خاک منجر به دستیابی به اطلاعات ارزشمندی می‌گردد که در اتخاذ تصمیمات مدیریتی برای فعالیت‌های 
کشاورزی موثر است. اين مطالعه با هدف بررسی تغییرات مکانی قابلیت هدایت الکتریکی (:6) به عنوان بيانگر وضعیت شوری خاک و کرین‌آلی خاک 
(906) در دو لایه سطحی (صفر تا۳۰ سانتی‌متر) و زیرسطحی (۶۰-۳۰ سانتی‌تر) با استفاده از چهار الگوریتم یادگیری ماشین جنگل تصلدفی (*8۳ 
درخت تصمیم (1211» رگرسیون بردار پشتیبان (۷۴) و شبکه عصبی مصنوعی ((۸) در دشت قزوین انجام شد. بدین منظور, از ۲۷۸ خاکرخ 
مطالعاتی نمونه‌برداری گردید و پس از عبور از الک دو میلی‌متر مقادیر ۳ و 5006 آنها اندازه گیری گردید. متفیرهای کمکی مستخرج از مدل رقومی 
ارتفاع و ماهواره لندست ۸ شامل ارتفاع (16۷۵1100) تابش پخشیدگی (1111056)» شاخص همواری کف دره با درجه تفکیک بالا 1۷۳۷2۳ شاخص 
تفاوت پوشش گیاهی نرمال شده (019۷1)» شاخص خیسی ساگا (91۷1) و شاخص جهت باد (1۷15) به‌عنوان نمایندگان عوامل خاکسازی مورد استفاده 
قرار گرفتند. نتایج نشان داد که مدل 1*۳ برای پیش‌بینی ۳ خاک در لایه سطحی با 122 برابر با ۰/۷۴ و 1۷۳ ۰/۳۶ و 012۷5۳ ۰۷/«وهمچنین 
پیش‌بینی کربن‌آلی در هر دو لایه سطحی و زیر سطحی با 157 به‌ترتیب برابر با ۰/۹۰ و ۰/۸۰ و به دنبال آن مدل 15 برای پیش‌بینی شوری خاک در 
لایه زیرسطحی با 122 ۰/۷۷ و /1۷]5۳ ۰/۹ و 0120۷5۳ ۰/۱۷ نسبت به سایر مدل‌ها دارای صحت بالاتری بودند. همچنین» متغیرهای توپوگرافی 
شامل ارتفاع شدت تابش پخشیدگی و شاخص همواری کف دره با درجه تفکیک بالا بیشترین ارتباط را با تغیبرات شوری و کربن‌آلی خاک در هر دولابه 
سطحی و زیرسطحی داشتند. به طور کلی مدل‌های "15۳ و 171۲ به همراه متغیرهای توپوگرافی توانستند تغیبرات شوری و کرین آلی سطحی و زیرسطحی 
خاک را با صحت قابل قبول در منطقه مورد مطالعه ارائه نمایند؛ که نقشه‌های تهیه شده می‌توانند برای اعمال تصمیم‌های مدیریتی لازم در مورد خاک‌های 
منطقه مورد استفاده قرار گيرند. 


واژه‌های کلیدی: تغییرات سطحی و عمقی» متغیرهای محیطی» مدلسازی رقومی» نقشه‌برداری ویژگی‌های خاک 


۱- دانشجوی دکتری پیدايش و رده بندی و ارزیابی خاک گروه علوم و مهندسی خاک واحد علوم و تحقیقات دانشگاه آزاد اسلامی» تهران» ایران 
۲- استاد گروه علوم و مهندسی خاک دانشکده کشاورزی, دانشکدگان کشاورزی و منابع طبیعیء دانشگاه تهران» کرج» ایران 
(#- نویسنده مسئول: ۱.۵6.1۳ ۵ 15271020 تلتقصط) 
۳- استاد گروه علوم و مهندسی خاکواحد علوم و تحقیقات. دانشگاه آزاد اسلامی » تهران» ایران. 
۴- استادیار پژوهشی موسسه تحقیقات خاک و آب. سازمان تحقیقات. آموزش و ترویج کشاورزی کرج» ایران 
۵- استادیار گروه علوم و مهندسی خاک واحد علوم و تحقیقات دانشگاه آزاد اسلامی» تهران» ایران 
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۸ نشریه آب و خاک. جلد ۰۳۷ شماره ۲. خرداد - تیر ۱۴۰۳ 


مقدمه 

بیش از یک سوم اراضی جهان تحت تاثیر شوری قرار دارند که 
رشد و عملکرد محصولات را تهدید می‌کند و مانع از توسعه کشاورزی 
پایدار می‌شوند (2020 ,. 6 220). شوری خاک یکی از بزرگترین 
چالش‌ها در مناطق خشک و نیمه‌خشک است که به شدت بر تولیدات 
کشاورزی تاثیر می‌گذارد (2019 .۵7 6۶ ۲۱211102). 

کربن آلی خاک (506) نقش مهمی در ظرفیت نگه‌داری آب» 
تبادل کاتیون-آنیون و در دسترس بودن مواد غذایی دارد که باعث رشد 
گیاه می‌شود. بنابراین» کاهش محتوای ۹06 در خاک های کشاورزی 
به طور کلی تهدیدی بزرگ برای پایداری کشت خاک درنظر گرفته 
می‌شود (2019 ,.۵1 6 217212). فرآیندهایی مانند زهکشی " و تیغه 
کشی " بسته به نوع خاک» سیلاب و شوری می‌توانند منجر به 
اکسیداسیون مستقیم 506 شوند (2018 ,.ا۵ ۶ اکنناوصام۳) 
بنابراین» بررسی خاک‌ها از نظر شوری و کربن‌آلی و توزیع مکانی آن 
برای حمایت از تصمیم گیرندگان در برنامه‌ریزی توسعه کشاورزی برای 
کاهش چالش‌های مربوط به امنیت غذایی در اراضی متاثر از نمک از 
اهمیت بالایی برخوردار است. اين در حالی است که گزارش وضعیت 
منابع خاک جهان (2015 ,1۳۳5۰ ۵00 ۳۸0۵) اشاره می‌کند که 
اطلاعات درمورد میزان و ویژگی‌های خاک‌های آسیب دیده از نمک 
بسیار پراکنده است. 

در سال‌های اخیر تعداد زیادی از مطالعات شوری خاک ( ,00۳۷۷10 
21 ,۸ 61 ۷۷۵8۵ :2021 ,.۵1 6 ۷۵ :2020) و محتوای کربن 
آلی ۱ 19 .۵۲ ۸ 00022060 :2019 ره 2 7000 
90 ,۵ ۵ 206طتقآ) مبتنی بر سنجش از دور انجام شده 
است. حسنی و همکاران (2020 ,./۵ 6 [۲12:520) با استفاده از مدل- 
های یادگیری ماشین ۷۲](۴) و محموعه‌ای از داده‌های اقلیمی. 
توپوگرافی, خاک و سنجش از دور 185(۵) به پیش‌بینی مکانی خاک‌های 
شور و سدیمی در عمقق خاک و با دوره‌های زمانی مختلف در مقیاس 
جهانی پرداختند. تقدسی و همکاران (2018 ,اه ۶ 20051ع12) با 
استفاده از تصاویر سنتینل ۱ و با کمک مدل رگرسیون بردار پشتیبان ۶ 
(۶۷) به پايش شوری خاک با 182 برابر با ۰/۹۷ و 115 برابر با 
۵ در اراضی با خطرشوری شدید پرداختند. حمزه‌پور و همکاران 


(2019 ,.اه 21 ام۳12۲۱26۲) مطالعه‌ای با هدف ارزیابی برهمکنش 
و تأثیر نسبی متغیرهای محیطی انتخاب شده بر تغییرات مکانی 506 
انجام دادند. نتایج نشان داد تعامل قابل توجهی بین متغیرهای ]2۷ 
٩2۸0 5‏ و تغییرات ۹06 وجود دارد. عبدالعماد و همکاران 
(2019 .۵1 ۶1 ۵۸0-۳[02000) به بررسی گسترش شهرسازی و تأثیر 
آن بر پتانسیل کشاورزی و ذخایرکربن آلی ب استفاده از دده‌های ماهواره 
ای چندزمانی و تکنیک‌های سنجش از دور در دلتای نیل پرداختند. 
نتایج حاکی از آن بود که گسترش شهری منجر به از دست رفتن ذخایر 
0 و تفییر خاک‌های حاصلخیز کشاورزی به خاک‌های شور و با 
زهکشی ضعیف شده است. 

آگاهی از توزیع مکانی شوری و ۹06 برای دستیابی به اطلاعات 
در اتخاذ تصمیمات مدیریتی برای کشاورزی پایدار موثر است. در 
اکثرمطالعات انجام شده به روش کمی» سنسورهای سنجش از دور برای 
برآورد شوری خاک در سطح به‌ کار گرفته شده‌انده این در حالی است 
که شوری خاک زیر سطحی بر رشد گیاه تاثیر گذار است زیرا مانع از 
جذب رطوبت و مواد مغذی توسط ریشه می‌شود ( 200 1۲10201 
1 ,110/271). العزیز و همکاران (2018 ,.۵ 2 ۵2۱2 ۳1طه) 
اعلام کردند که اراضی دارای ۳ بین ۲/۲۵ تا ۵ دسی‌زیمنس بر 
سانتی‌متر از ۱۷ درصد در سال ۱۹۹۱ به ۳۷ درصد در سال ۲۰۰۶ 
افزايش یافت. سوالیم وهمکاران (2019 ,./0 61 )٩۱/1162۳0‏ نقشه‌های 
رطوبت و شوری خاک را در عمق ۳۰ و ۶۰ سانتی‌متری خاک تهیه 
کردند و نشان دادنده که ضریب همبستگی بین شوری و رطوبت لایه 
سطحی و زیر سطحی به‌ترتیب ۰/۴۸ و ۰/۶۸ بود. آولیاکولف و همکاران 
(2020 ,.۵1 ۵1 ۵۷11۱/26010۷) با استفاده از داده‌های سنحش از دور 
و تکنیک آنالیز نقشه‌بردار زاویه طیفی (۹۸۵) با دقت کلی ۲۴/۷۴ 
درصد نقشه خاک‌های متاثر از نمک را در اعماق مختلف تهیه کردند. 
اردکانی و وحدتی (2017 ,۷۵502 0صه تصفآد0ته) تغییرات 506 
و شوری خاک را در بین سال‌های ۲۰۰۸ تا ۲۰۱۲ در استان یزد مورد 
بررسی قرار داده و نتیجه گرفتند که در طی این مدت محتوای کربن‌آلی 
اراضی کشاورزی افزایش يافته و مساحت اراضی شورکاهش يافته است. 

با این حال, مطالعات بسیار کمی به بررسی همزمان تغییرات 
سطحی و زیر سطحی تغییرات مکانی قابلیت هدایت الکتریکی ۳0۳ 
به عنوان بیانگر وضعیت شوری خاک و محتوای کربن‌آلی خاک در 


1- ٩011 022016 20 

ععفقصنه0 -2 

عصلانل -3 

عصتصنجم 1 مصتطم۱۷]2 -4 
8 1601016 -5 

مصتط)2] ۱۷ ۷۵۵۲۵۲ 9۱۵۵۲۲ -6 
۵۵۵۲ 20216 50601۲21 -7 
۲ ۳۱۵۵1621 -8 


رضایی و همکاران. نقشه‌برداری رقومی تغییرات سطحی و زیرسطحی کربن آلی و شوری خاک در بخشی از اراضی دشت قزوین ... ۳۱۹ 


اراضی کشاورزی پرداخته اند. اکثر مطالعاتی مبتنی بر تغییرات سطحی 
ویژگی‌های خاک از جمله کربن آلی و شوری در سطح معطوف بوده 
(2017 .۵1 ۵ نفک :2020 .۵ 61 تلهة1) و در کمتر مطالعاتی 
به بررسی همزمان تغییرات سطحی و زیر سطحی آنها پرداختند. با توجه 
به اينکه جمله شوری یک ویژگی دینامیک خاک است و تغییرات آن در 
دوره‌های مختلف زمانی از سطح به عمق متغیر است. با توجه به تاثیر 
متقابل شوری و کربن‌آلی برروی سایر ویژگی‌های خاک و همچنین 
سایر کارکردهای اصلی خاک اهمیت این بررسی فزونی می‌یابد. از دیگر 
سوء با توجه به اینکه بخش عمده منطقه مورد مطالعه تحت کشت 
گیاهان یکساله می‌باشد که نفوذ و توسعه ريشه آن‌ها تا اعماق 
زیرسطحی نیز گسترش می‌یابنده درک و اطلاع از تغییرات مقلیر 
شوری و کربن آلی خاک می‌تواند بر مدیریت تغذیه. حفظ حاصلخیزی 
و افزايش عملکرد محصولات زراعی موثر باشد. بنابراین جهت پوشش 
بیشتر اطلاعات در این زمینه مطالعه حاضر با اهداف ۱) تهیه نقشه 
پیش‌بینی مکانی شوری و کربن آلی خاک در دو لایه سطحی (۲۰-۰ 
سانتی‌متری) و زیر سطحی (۶۰-۲۰ سانتی‌متری) با استفاده از چهار 
الگوریتم یادگیری ماشین جنگل تصادفی (8۳» درخت تصمیم 
(011) رگرسیون بردار پشتیبان (5۷8) و شبکه عصبی مصنوعی 
[(۸۱) و ۲) شناسایی مهمترین متفیرهای کمکی محیطی تاثیرگذار 
بر روی پیش بینی مکانی شوری و کربن آلی خاک در دو عمق سطحی 
و زیرسطحی مورد نظر در بخشی از اراضی دشت قزوین که تحت کشت 
و برداشت دائمی بوده و از توان تولید قابل توجهی در تولید محصولات 
کشاورزی برخوردار است طراحی گردید. 


مواد و روش‌ها 

منطقه مورد بررسی شامل بخشی از اراضی کشاورزی و عرصه‌های 
طبیعی استان‌های البرز و قزوین. حدفاصل شهرستان‌های نظراباد و 
جغرافیایی ۵۴ ۲۵۶ تا ۵۴ ۲۶۶ شمالی و طول جغرافیایی ۱۵ ۵۰۴ تا 
۰ شرقی واقع شده است (شکل ۱-لف» ۱-ب و ۱-ج). حداقل 
و حداکثر ارتفاع در این منطقه به ترتیب ۱۱۴۰ متر و ۱۷۴۷ متر نسبت 
به دریای آزاد است. این منطقه دارای اقلیم خشک و نیمه‌خشک 
براساس روش کوپن است. واحدهای فیزیوگرافی منطقه مورد مطالعه 
شامل تیه و فلات در قسمت شمال» دشت دامنه‌ای قسمت مرکز و 
اراضی پست در جنوب می‌باشد. مواد مادری تشکیل دهنده خاک‌های 


منطقه بر اساس نقشه زمین‌شناسی (۱:۱۰۰۰۰۰) سازمان زمین‌شناسی 
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کشور شامل رسوبات آبرفتی و کوهرفتی مربوط به دوره کواترنر بیشتر 
متمرکز در بخش‌های شمالی و مرکزی منطقه و کفه‌های گلی و نمکی 
عمدتا در بخش‌های جنوبی منطقه می‌باشد. خاک‌های منطقه در سه 
رژیم رطوبتی زریک خشک» اریدیک ضعیف و آکوئیک و حرارتی 
ترمیک واقع شده‌اند (2022 ,اه 6 ت۳۵2) و از نظر رده‌بندی در 
سه رده اریدی‌سولز, انتی‌سولز و اینسپتی‌سولز قرار دارند ( ٩0۳۷6‏ 5011 
7 ,/2اه). 


روند اجرای تحقیق 

شکل ۲ روندنمای اجرای تحقية را در۴ مرحله نشان می‌دهد: ۰۱ 
حفر نقاط نمونه‌برداری و اندازه‌گیری خصوصیات فیزیکی و شیمیایی 
خاک و تهیه متغیرهای کمکی, ۲. مدل‌سازی پراکنش مکانی ۳ و 
درصد 906 در لایه سطحی و زیر سطحی توسط مدل‌های یادگیری 
تعیین آهمیت نسبی متفیرهای محیطی و ۴. تهیه نقشه پیش‌بینی توزیع 
مکانی سطحی و زیر سطحی 9 0 در منطقه مورد مطالعه 


نمونه‌برداری و تجزیه و تحلیل آزمایشگاهی 

مظا اف میتی ۵ یاس رای شین تا احان رن 
یک دوره ۶ سال) انجام شد و براساس روش طبقه‌بندی تصادفی در 
نهایت ۲۷۸ خاکرخ مشاهداتی مورد مطالعه واقع گردید (شکل ۱-چ). 
نقاط مشخص گردیده تا عمق ۱۵۰ سانتی‌متری يا بالاتر از لایه 
محدودکننده حفر و براساس راهنمای تشریح صحرایی خاک‌ها تشریح 
و از کلیه افق‌های ژنتیکی قابل شناسایی نمونه‌برداری شد. پس از 
نمونه‌برداری و شماره‌گذاری» نمونه‌های خاک به آزمایشگاه منتقل شد. 
پس از هوا خشک شدن و عبور تمامی نمونه‌ها از الک دو میلی‌متر, 
مقدار قابلیت هدایت الکتریکی خاک (0ظ) در عصاره گل اشباع 
(1982 ,:50206) به‌عنوان بیانگر وضعیت شوری خاک و کری ن‌آلی 
خاک به روش واکلی بلک (1975 ,1201508) برای تمامی نمونه‌ها 
اندزهگیری شد. در نهایت میانگین وزنی 36 و 506 برای اعماق ۰- 
۰ ۶۰-۳۰ سانتی‌متری محاسبه شد. 


متغیرهای محیطی و انتخاب آن‌ها 

در نقشه‌برداری رقومی خاک توزیع مکانی کربن آلی و شوری خاک 
و روابط آنها با عوامل تشکیل‌دهنده خاک از جمله خواص خاک مواد 
مادری» اقلیم» موجودات زنده. پستی و بلندی و زمان براساس مدل 
مکانی اسکورپن که توسط مک برتنی و همکاران ( ۶ ۱۷۲۵۳۲۵06۷ 
3 ,.0) ارائه شده توضیح داده می‌شود. 
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شکل ۱- الف. نقشه استان‌های قزوین و البرز؛ ب و ج. منطقه مورد مطالعه و پراکنش مکانی نقاط نمونه‌برداری شده 
ععام‌صصوی آزمز ۵۶ حمتاتحاتصاعتل اهتادمری 0ص مه واه ]۵۱ مممصز ماهر 002 0 002210 م۸ -1 ما۲ 


متغیرهای توپوگرافی از مدل رقومی ارتفاع (آلوس پالسار ) با وضوح 
۵ متر در نرم افزار 075 5۸0۸ مستخرج و متغیرهای سنجش از 
دور از میانگین تصاویر مربوط به فصل‌های تابستان» پاییز و زمستان 
(ماه‌های شهرپور تا بیمن) طی سال‌های ۱۳۹۵ 2 ۱۳۰۰ ماهواره 
لندست ۸ (با شماره ردیف و گذر به ترتیب ۲۵ و ۱۶۵ با قدرت مکانی 
۰ متر) با استفاده از نرم افزار] ۳۲۷ نسخه ۵,۲,۵ تهیه شدند. ابتدا 
کلیه پیش پردازش‌های لازم از جمله تصحیحات رادیومتری و بهبود 


ص۲۳۸ ونم -1 


کنتراست انحام شد. در این تحقیق متغیرهای محیطی مورد استفاده 
پراساس مطالعات پیشین صورت پذیرفته توسط رحمانی و همکاران 
(2022 ,.۵ ۵۶ 2۳021)» موسوی و همکاران ( ,.۵1 ۶ 1۷0952۷1 
2 و خاموشی و همکاران (2017 ,۵ 6 تحاد0تحطک) در منطقه 
مورد مطالعه» برای پیش‌بینی مکانی ۳ و ۹06 انتخاب شدند (جدول 


0 
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۱ 5 تصاویر ماهواره لندست ۸ مدل رقومی ارتفاع 
۲ ممناه۱۵7ظ آمازع‌نط 
لارد ( ۳۰-۰ مساقشی مترم 27 ۶ )482‌صما ه 
لایه زیر سطحی (۶۰-۳۰ سانتی متر) توودت 
1 
صصعصصعتعممص ۳۵ هه 506 با طصیح ی پارامترهای تو و گرا: 
(حصه 0-30) عبرعا آنمو ورو۲" سنجش از دور پارامترهای ثوپو تراگی 
(جصت 30-60) م1 آتمو ماناه عصتعده5 عامحصع۳ 10002۲۵ 
۳ ههور 


سس متقیرهاین کمتکی تتفاي.شتاة 
سرا هتسخ 5۵160160 
اوری ۹۱0۰۵۵ ۱5 و 
مدل سازی مکانی گام ۱ 
متغیرهای کمکی + 506 + 170 1 500 
۱ ۱ جمع آوری داده 
۵8۵ 5۵0121 مدع د)دط 
وماطاهته۷ دتا ۵9۲۵ وب ۲ 
گام ۲ 
۳ ۳0 1 2 560 
دوه میج ب یی | اد اه | | یل عناق انتخاب مدل های 
۳ | ی از س || 9 یادگیری ماشین 
666۵ ۵0615 ]۷۲ 


مناسب ترین مدل پیش بینی مکانی ِ 5 گام ۳ 
61مص صمت)ین0عم اهتاهوو افعظ ماوقا کی 503 
۶ ععصهاتمجرحصز ۲۹۵۱۵۷۵ صحت سنح 
۳ ,۳ ]1۱۷ ,1*2 ععاطهنته۷ حهتاتسه ۷ 
نقشه توزیع مکانی 50 نقشه توزیع مکانی ۲6 گام ۴ 
ی : ان 2 ۰ 4 560 
۶ طمنباطا تون نادمه ۶ )اما تون آهنادم5 نهیه نقشه 
مقصط 5۵ ۱۱۱ ممناددوعم ود 


شکل ۲- روندنمای پیش‌بینی مکانی شوری و کربن‌آلی خاک 
۸ ۵۵۱۵و ممصای متصجعن 0صه واتصتامی ازمی ۵۴ )تهم ۲۱۵۳۲ -2 ما۲ 
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جدول ۱- متغیرهای کمکی مورد استفاده در پیش‌بینی ۳ و 50 
60 20 ۲ م)ز0ع۱ص نز 0عون وماه1 ۷۵۲ هتاتت۸ -1 ماطاو1 


نماینده عامل 5 1 
شرح نماد 
خاکساز ۳ 
۱/9۱ مر 
ار( 
توپوگرافی 
4 
1 تابش مستقیم انعکاس يافته از سطوح 
توپوگرافی "قلیم 4 ان 
۳ حم 1۲61166160 مرمنا ۲201 11۲۵6۵۲ ۳6 
286 
توپوگرافی-قلیم شاخص رطوبت توپوگرافی بر اساس حوضه آبریز ۳ 
1-0 اصلاح شده 
توپوگرافی"اقلیم مناطق در معرض یدبا توجه به جهت با 1۳ 
1-0 
توپوگرافی اندازه‌گیری همواری و پستی و 
1 امه همه آقتتاووم صع وعمصه[۳۲ 
موجودات-خاک رو 0 
تن وربور 


ژئومورفومتری مستخرج از مدل رقومی ارتفاع می‌باشد و همانطور که 
از نامش مشخص است. شبیه به «شاخص خیسی توپوگرافی» است. 
اما بر اساس یک حوضه آبریز اصالاح شده محاسبه می‌شود که به مقدار 
کمتری بیانگر میزان جریان می‌باشد. در نتیجه برای سلول‌های 
(پیکسل‌های) واقع در موقعیت‌های پست زمین‌نما با فاصله عمودی کم 
تا شبکه آبرهه یا کنال مجاوره رطوبت خاک را با پتانسیل واقی‌تری 
در مقایسه با شاخص استاندارد خیسی تویوگرافی محاسبه می‌نماید؛ که 
اما بل خاوت هی خاشصی خی اک و خاخضی ی یوگ اف 
ماژول مجزا تهیه می‌گردند. در رابطه روبه رو معادله ریاضی محاسبه 
شاخص خیسی ساگا (۹۷۷1) ارائه شده است ( .۵1 6 ۳06006۲ 
02 


۸ -- 
7 صه/ ( وت 
در اینجا: 5۷۷1 شاخص خیسی ساگ؛ 150۸ < حوضه آبخیز 
اصلاح شده و 2 20 تانژانت زاویه شیب بر حسب رادیان می‌باشد. 


مدل‌سازی مکانی 


چهار مدل جنگل تصادفی, ماشین بردار پشتیبان, درخت تصمیم و 
شبکه عصبی مصنوعی با استفاده از بسته‌های تخصصی ۲2۸00۱ 
امعت۳ طاحاصمکآ دای ]12۲ تعصط و 027 و کد نویسی در نرم‌افزار 
0 نسخه ۴,۱,۱ برای پیش‌بینی توزیع مکانی ۳ و ۹06 در 


دو عمق سطحی و زیر سطحی مورد استفاده قرار گرفت. به طور کلی 


«مناعص۳۲ -1 
۵ 2712016 ۷ -2 


۱ 


عم۲2۱۵0 هتاتسش۸ 0۵۰۵ 
ارتفاع 
۱۳ 


شدت تابنش پخشیده 
ماه ۱۱/۱۱۱۱۱۱۱۵ 


۳ مدل رقومی ارتفاع 
شاخص خیسی ساگا 
ات ۲ (1]2م1([ 
۷۷۵0699 ۸۸ م00 
/۱0 ۲7۷ 
شاخص همواری کف دره با درجه تفکیک بالا 
متام ما۷ ۵۲ 1۵06 ممتاتاا0عم؟- ات۷۲ 
وععص۳[2 
شاخص تفاوت پوشش گیاهی نرمال شده سنجش از دور 


۴ (686)۵010 ۲۷ ۵۲۵86۵5]]ز12 ۱۱۵۱9۵۱1260 عط 6 مامصصم؟ 


مجموعه داده‌ها 

به صورت تصادفی ۸۰ درصد برای آموزش و ۲۰ درصد برای 
ارزیابی مدل‌ها تقسیم گردیدند. درجه اهمیت هریک از متغیرهای 
محیطی در مدل‌های یادگیری ماشین مورد استفاده طبق تابع۲ اهمیت 
متغیر ( ورام در بسته 02761 انحام گردید. همچنین جهت حصول 
پایداری نتیج در برازاش مدل‌ها از رویکرد ۱۰ گام ۱۰ تکرار استفاده 
گردید. نحوه دقیق عملکرد مدل‌ها درپیش‌بینی ویژگی‌های ذکر شده به 
شرح ذیل ارائه شده است: 


جنگل تصادفی (18) 

جنگل تصادفی که توسط بریمن (2001 ,1361۳0۵0) توسعه یافته 
استه یک روش یادگیری ماشین است که به طور گسترده بای طبقه- 
بندی و رگرسیون استفاده می‌شود. یک درخت تصمیم منفرد به دلیل 
واربانس و سوگیری بالای آن طبقه‌بندی و رگرسیون ضعیفی را ارائه 
می‌دهد. در حالی که ۳ تمایل به تولید مدل‌های قوی دارد زیرا می - 
تواند مشکلات یاد شده را با استفاده از مجموعه‌ای از درختان کاهش 
دهد. 1۳ هزاران درخت دوتایی تصادفی را برای ایحاد جنگل تصادفی 
تولید می کند. هر درخت براساس یک نمونه بوت استرپ. با استفاده از 
روش درختان طبقه‌بندی و رگرسیون (۸8۲)) با زیرمجموعه 
تصادفی از متغیرهای انتخاب شده در هر گره رشد می‌کند. برای هر 
درختی که بر روی یک نمونه بوت اسرپ رشد می‌کنده نرخ خطای خارج 
از کیسه (0013) با استفاده از مشاهدات خارج از بوت استرپ محاسبه 


68 1۵21۲095100 مه ط0تا2ع تمعن -3 
2 0۴ 0۷9 -4 


رضایبی و همکاران. نقشه‌برداری رقومی تغییرات سطحی و زیرسطحی کربن آلی و شوری خاک در بخشی از اراضی دشت قزوین ... ۳۳۳ 


می‌شود. برای اجرای مدل 1*۳ کاربر باید دو پارامتر پیش‌بینی» تعداد 
درختان در جنگل (0066) و تعداد متغیرهای قابل استفاده در انتهای هر 
گره (1۳۷) را بهینه کند تا خطای 0053 را به حداقل برساند و عملکرد 
مدل خوبی را به دست آورد. 


درخت تصمیم (1۲) 

مدل درختی تصمیم‌گیری طبقه‌بندی» یک روش ناپارامتری است 
که قادر به پیش‌بینی متغیرهای کمّی یا متفیرهای طبقه‌بندی شده بر 
اساس مجموعه‌ای از متغیرهای پیش بینی کننده کمّی و کیفی است. 
رشد درخت با پیوند و انشعاب و مرحله دوم. مرحله توقف و هرس 
می‌باشد. هدف از مرحله دوم به حداقل رساندن خطای پیش بینی است 
فرآیند طبقه‌بندی یک ویژگی از ريشه درخت شروع می‌شود داده‌های 
آزمون ارزیابی شده 9 شاخه مناسب برای نتبحه انتخاب می‌شود 9 تا 
افاتیکت یک وگ طای که ای قرا شم اکامهم انش با طاه شخ 
برگ قفا شود که آن وف سل به کلاسی فامگتارین شاه بط 
برگ است. روش درخت تصمیم برای طبقه‌بندی کاربری اراضی دارای 
مزایای قابل توجه است زیرا بین ویژگی‌ها و طبقات انعطاف‌پذیری و 
توانایی مدیریت روابط غیر خطی در این مدل وجود دار بنابراین دقت 
طبقه‌بندی را افزايش می‌دهد (2018 .اه 61 تصحصطدل). 


شبعه عصبی مصنوعی (۸(۷) 

شبکه‌های عصبی مصنوعی گروهی از الگوریتم‌های غیرخطی و 
مبتنی بر داده‌های عمومی هستند که پس از قرار گرفتن در معرض 
تغییرات ناگهانی پایدار می‌مانند و می‌توانند در محیطهای غنی از داده, 
آموزش داده و تولید شوند ( ,۵1 6 نمهزیم0ظ-سممنمنل‌مهکوظ 
0 ساختارهای شبکه عصبی مصنوعی ((۵) از مجموعه‌های 
به هم پیوسته نورون‌ها تشکیل شده‌اند که همبستگی‌های غیرخطی بین 
هر متغیر را تخمین می‌زنند. متغیرهای محیطی به عنوان نورون‌های 
ورودی ابتدا به یک پا چند لاه از نورون های پنهان و سپس به نورون 
های خروجی که متغیر هدف خاک هستند متصل می‌شوند ( ۲۲۵0۳۵ 
6 ر.6 61). 


ماشین بردار ی تید ۳ (۷) 

ماشین بردار پشتیبانی یکی از روش‌های یادگیری ماشین معرفی 
شده توسط بوسر و همکاران (1992 .۸ 2 05:6() است که براساس 
نظریه یادگیری آماری توسعه داده شده است. مدل‌های ماشین بردار 


طصتصصهاه 1۳[ اصعلم006111 -1 
۲ 500۵۲۵ 29ع]۱۷ ]500 -2 


کننده ماشین بردار پشتیبان و ب: مدل رگرسیون بردار پشتیبان. برای 
حل مسائل طبقه‌بندی داده‌ها در کلاس‌های مختلف از مدل ماشین 
بردار پشتیبان و برای حل مسائل پیش‌بینی از مدل 5۷15 استفاده می 
شود. منظور از رگرسیون به دست آوردن ابر صفحه‌ای است که با داده- 
های داده شده برازش می‌کند. فاصله از هر نقطه در اين ابر صفحه 
نشان‌دهنده خطای آن نقطه خاص است (2020 ,4 6 نصدحا2ط). 
صحت‌سنجی 

فرایند اعتبارسنجی مدل‌های تولید شده با استفاده از ۸۲۰ از کل 
داده‌هاء انجام گرفت. برای این منظور مقادیر برآورد شده برای هر نمونه 
از زیرمجموعه‌ای مستقل تعیین شد. صحت مدل‌ها با استفاده از شاخص 
های آماری ضریب تعیین ۱ (*18» میانگین ريشه مربعات خطا"(801515]) 
و میانگین ريشه مربعات خطای نرمال شده" (0180۷1518) برآورد شد. 
رتبه بندی مقادیر ضریب تبیین (معادله ۱) برای پیش‌بینی ویژگی‌های 
خاک عبارت است از: خیلی خوب (18<0.81» خوب (-0.61 < 1 
1 متوسط (0.410.6 - 8) و ضعیف (482>0.4 (2009 
,۵۷ وه ۴8055۵1). ۷8۳ یک معیار پرکاربرد برای 
تخمین تفاوت بین مقادیر پیش‌بینی شده توسط یک مدل و مقادیر 
اندازه گیری شده است. همچنین :190۷15۳ مقیاسی از انداره خطا را برآورد 
می‌کند که به مقادیر پرت حساس است و توسط معادله ۲ محاسبه می 
شود. همچنین مقادیر ۱۷۲5۳۲ (معادله ۲) کمتر از ۰/۱ نشان دهنده 
پیش‌بینی عالی» بین ۰/۱ و ۰/۲ خوب. بین ۰/۲ تا ۰/۳ متوسط و بیشتر 
از ۰/۲ ضعیف خواهد بود (2006 ,.اه 6 ۷۷۵112 ). 
(-۷۱-۲(0۱) ,۲۱۳ 


معادله ۱ < 2[ 
۴۱-2) ,9 0۱-۰2 
اه ۱ 
۷ 


در اینجا :7 داده اندازه‌گیری شده و7 مقدار پیش‌بینی شده است. 
همچنین 7 میانگین داده‌های اندازه‌گیری شده و میانگین داده‌های 
پیش‌بینی شده است. همچنین اه تعداد داده‌ها است ,.4 6 هاعاظ۳ 
(2020). 


معادله ۳ 0« 


(22 


۳۴(۷]9] < 


۲ 0۱۷۵7 2ع]۱۷ ]00 ۱0۲۴۱۵11260 -3 


۴ نشربه آب و خاک. جلد ۳۷. شماره ۲. خرداد - تیر ۱۴۰۲ 


جدول ۲- توصیف آماری شوری و کربن‌آلی در دو لایه سطحی (۲۰-۰ سانتی‌متر) و زیر سطحی (۰۰-۳۰ سانتی‌متر) 
12۵5 رن 630-60 عهاو-طایاز صرح رصن 0-30) هساو ۳۵ صذ مدای »تصهعه 0ص انصناهه ۵ ی‌تامناهاه ۱(۵50۵/0)۱۷۵ -2 12۳016 


انحراف معیار میانگین حداکثر واحد حداقل ویژگی‌های خاک 
06۷126100 0270 ۱/4 ]۱۷ نصا ص۱۷ ناه روم ام 
لایه سطحی (۲۰-۰ سانتی‌متر) 
ره 0-30) تعاع1 501 100 
102 05 19111 دسی‌زیمنس بر متر 014 هذایت الکتریکی (۱۳۳ 
5/1 (0ظ) تالجم لحعت00ع۱ظ 
و0 084 33 99 011 یی ی خاک ۱۳۰ 
9 (800) طمطاعدن منصقع:0 5011 
لایه زیر سطحی (۶۰-۳۰ سانتی‌متر) 
(ت 30-60) م1 لر0ووانن 
1443 555 1085 دسی زیمنس پر متر 008 ری ۱۳۳ 
/ (۱0ظ) تالجم لحعت00ع۲۱ 
029 04 209 0 002 یآ ۱۳۳۰۳۵ 
10 (00) ممطایدن منصهع:0 011 
نتایج ۶ و ۲۹/۴۷ و متغیر 1111156 ۲۸/۱۵ و ۲۷/۶۷ درصد از تغییرات 
۲ ‌ 0 (شکل ۲ج و ۲-د) ر در لایه سطحی و زیر سطحی را نشان 
توصیف اماری 


در جدول ۲ ارائه شده است. مقادیر حداقل و حداکثر 76 در هر دولایه 
سطحی و زیر سطحی نشان داد که منطقه مورد مطالعه دارای خاک‌های 
غیر شور تا بسیار شور تا عمق ۶۰ سانتی‌متری است. همچنین, ميانگین 
دسی زیمنس بر متر از لایه سطحی به لایه زیر سطحی نشان داد. اين 
یافته‌ها با نتایج فتحی‌زاده و همکاران (2020 ,.اه 6 12۵060ظ۳۵۱) 
منطبق است. میانگین 906 از سطح به عمق به ترتیب برابر با ۰/۸۴ 
و ۰/۴ درصد است. این در حالی است که بیشترین مقدار ٩06‏ در لایه 
سطحی و برابر با ۳/۲ درصد مشاهده شد. براساس طبقه‌بندی ویلدینگ 
(1984 ,ع۷۷11010) کلاس تغییرپذیری ۳ و 806 در هر دو لایه 
سطحی و زیر سطحی در کلاس حداکثر (بیش از ۳۵( قرار دارد. 
تغییرپذیری زیاد ۹06 و 6 خاک نشان داد که خاک‌های منطقه مورد 


مطالعه تحت تاثیر مدیریت» متفیرهای محیطی و تخریب قرار دارند 
(2020 ,.۵1 ۶۶ طع0220نامحصطه۱۷). 


متغیرهای کمکی و اهمیت آن‌ها 

اهمیت متغیرهای کمکی استفاده شده برای پیش‌بینی ۵ 9۵ 
ی هت ی 
متر) در شکل ۲ نشان داده شده است. متغیر ۳16۷۵۸:02 و 1:56 
به عنوان پیش‌بینی کننده‌های مهم ۳ و۹0 در هر دو لایه سطحی 
و زیر سطحی تعیین شدند. متغیر ارتفاع ۵۰ و ۳۵ درصد و متغیر 
عمسلز< ۱۹/۷ و ۲۵/۱ درصد از تغییرات ۳ (شکل ۲-لف و ۲-ب) 
را در لایه سطحی و زیر سطحی تبیین نمودند. متغیر ۳۱6۷2100 


دادند. همچنین ۷۲۳۷۳۴۳ با ۱۹/۷ و ۲۰/۶ درصد اهمیت نسبی برای 
پیش‌بینی 6 در لایه سطحی و زیر سطحی, سومین متغیر کمکی میم 
بود. به طوری که مقادیر بالای ۳ در منطقه مورد مطالعه با مقادیر 
بالای 1۷۲۳۷2۳۲ منطبق بود. 

متغیرهای توپوگرافی به عنوان تاثیرگزارترین متغیرهای پیش‌بینی 
کننده برای پیش‌بینی هر دو متغیر ۳ و ۹06 در دو عمق سطحی و 
زیر سطحی مورد مطالعه شناسایی شدند. به طورکلی توپوگرافی یکی از 
مهمترین عوامل تشکیل خاک است و متغیرهای توپوگرافی مبتنی بر 
1۱<(۷ معمو لا به عنوان پیش‌بینی کننده‌های کلیدی برای رقومی کردن 
نقشه‌های خاک استفاده می‌شوند. پستی و بلندی‌ها با کنترل جریان آب 
و رسوبگذاری به نوبه خود بر توسعه خاک و توزیع مکانی خواص خاک 
تاثیر می‌گذارند (2020 ,.» 61 20۷). مطابق با یافته‌های تقی‌زاده و 
همکاران (2021 ,.01 ۶ 22۳01272060 ]) متغیرهای ۲۳۷۲ و 1۷۲۷۲۳۲ 
مهمترین‌متغیرهای کمکی در پش‌بینی ۳ خاک سطحی و زیر 
سطحی شناسایی شدند (شکل ۴). ژو و همکاران (2020 ,.41 6 20) 
نیز در مطالعه خود اعلام کردند که متغیرهای توپوگرافی دارای 
همبستگی بیشتری نسبت به سایر متغیرهای کمکی مورد استفاده با 
محتوای ۹06 و نیتروژن خاک داشتند. ژو و همکاران ( .۵1 ۶ 20 
9 نیز تاکید کردندکه متغیرهای توپوگرافی به عنوان مهم‌ترین 
متغیرهای پیش‌بینی کننده ۹06 در یک منطقه در یک بازه زمانی 
بیست ساله شناسایی شدند. در همین راستا محمودزاده و همکاران 
(2020 .۵1 ۶۶ ۷]2۵0۴0002206) نیز نتایج مشابهی را در پیش‌بینی 
5900 در غرب ایران گزارش نمودند. 
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شکل ۲- اهمیت نسبی متغیرهای کمکی برای پیش‌بینی الف: 86 لایه سطحی (۳۰-۰ سانتی‌متر) براساس مدل 18 ب: 80 لایه زیرسطحی (۲۰- 
۰ سانتی‌متر) براساس مدل 18۳ ج: »90 لایه سطحی (۳۰-۰ سانتی‌متر) براساس مدل 8۲ د: »50 لایه زیر سطحی (۱۰-۳۰ سانتی‌متر) 
براساس مدل 13۲7 

۲ 06 ۵ 0موو رصن 0-30) رها م۵ مظ ۵ ۱ بو 0ص بم مماصاهزه۱ هتاتده اه مصماته0 رازه۲۵ م1 -3 ما۲ 


۲ 06 ۵ 860و رصن 0-30) مرها (ز50 نم ۳6 ]۵ 9۹0۵۵ :۰ ما۵0 1۸۲ 6 جدن 0عقوها رصم 30-60) ۱/۵ ]تمد ماو 6 ۵۲ ۳۲ :۵ و1۵۵0 
۵06 1۵۳ ۳6 0 60عو رصصن 30-60) مرها ازم5 دبای 6 ۵۶ ۵۵ :02 و1۵۵0 


ارزیابی عملکرد مدل‌های یادگیری ماشین 

عملکرد چهار مدل 1۳ ۸۵ 1۲۳ و 5۷۴ با استفاده از 
متغیرهای کمکی منتخب برای پیش‌بینی ۳ و ۹0 در دو لایه 
سطحی و زیر سطحی (۲۰-۰ و ۶۰-۲۰ سانتی‌متر) با استفاده از 
شاخص‌های ارزیابی *12 ۳8۳ و 2۳۷۲5۳ در جدول ۲ ارائه شده 
است. نتایج پیش‌بینی 0 در لایه سطحی نشان داد که مدل 5۳ و 
۲ دارای بیشترین میزان صحت با مقادیر *12 به ترتیب برابر با ۰/۷۴ 
9 ۰/۷۲ و ۷۲ به ترئیب برابر با ۰/۳۶ و ۰/۸۲ دسی‌زیمنس بر متر 
می‌باشند. این در حالی است که در لایه زیر سطحی مدل 1717 با مقادیر 
1و برابر با ۰/۷۷ و ۰/۹ دسی‌زیمنس بر متر بیشترین صحت 


را برای پیش‌بینی ۳ داشت. به‌طور کلی هر دو مدل 3۳ و 117 با 
۳ بیش از ۰/۷ در هر دو لایه سطحی و زیر سطحی دارای صحت 
برای لایه‌های سطحی و زیر سطحی دارای عملکرد خوب برای پیش- 
بینی ۳ بود و مدل ٩۷۴‏ عملکرد متوسط داشت. همچنین نتایج 
پیش‌بینی ۹0 خاک در اعماق مورد مطالعه بيانگر این است که مدل 
۳ در پیش‌بینی کربن‌آلی خاک در هر دو لایه سطحی و زیر سطحی 
دارای بیشترین صحت با مقادیر ۳ برابر با ۰-۱۸۹ و ۰/۹ و 5۳( 
برابر با ۰/۱ ۰/۰۳ بود. با توجه به راسل و مکبرتنی ( 2040 561و 
9 ,/۱۷16۲۵106۷) سایر مدل‌ها با 187 کمتر از ۰/۴ عملکرد ضعیفی 
در نیشن نیت کرینآلی خاک داشتند 
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شکل ع- توزیع مکانی متغیرهای کمکی انتخاب شده الف:ار تفاع» ب: تابش پخشیده. ج: شاخص همواری کف دره با درجه تفکیک بالاء د: شاخص 


تفاوت پون سئیر ۱ گیاهی نرمال نیده» ۵: شاخص خیس ساگاء خ: شاخص اثر باد 
۶ ۲۱/6 1۸۵60۱0۵۸10- ۱۷۲۵۱6۱ ۰ مصمتاهاه‌عصا معسلزظ :0 رمتاه۲۱6۲ :۵ مماماهتیه۲ وهتاتسه ۵۱6060 ۱۴ متا ماتصاعنل امتادمرک -4 مساع۲1 
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مقایسه میانگین انحراف نسبی بین مقادیر مشاهده شده و پیش‌بینی 
شده برحسب درصد (ظ28[۷15) برای ویژگی‌های 130 و 806 نشان 
می‌دهد که مقادیر ۱11617 محاسبه شده توسط مدل 5۳ نسبت به 
فتایر فذال ها کش انشا براین سل ظخاک: در لاه سا 
و ۹06 در لایه زیر سطحی با ۱۷15 برابر با ۰/۰۷ با درجه دقت 
عالی و پیش‌بینی توزیع 12 در لایه زیر سطحی و 506 در لایه 
سطحی با درجه دقت خوب انجام شده است (2006 ,.۸1 6 21120 ۷۷). 
هنگل و همکاران (2004 ,.۸1 2 ۳1۵:816) بیان کردند که ۲۷۲5۳۲ 
بیش از ۰/۷ دقت پیش‌بینی غیرقابل قبولی دارد. بنابراین 3 دو 
محدودیت اصلی دارد که نمی‌تواند سوگیری پیش‌بینی‌های ضرایب و 
همچنین قابلیت اعتماد مدل رگرسیون را تعیین. به‌عنوان مثال با وجود 
چند نقطه تاثیر گذار یا پرت» **1 ممکن است برای مدلی با مقادیر خطای 
کوچک (مدل‌های مناسب) کم و یا برای مدلی با مقادیر خطای بزرگ 
(مدل‌های برازش ضعیف) زیاد محاسبه شود ( 204 صهن0ع121 
0 ,۷11116). ماشالابا و همکاران (2020 ,.ه ۶ 2+021202) 
گزارش کردند که صحت پیش‌بینی مدل 1 در عمق زیر سطحی با 
برابر با ۰۱۶/ بیشتر از عمق سطحی برای پیش‌بینی توزیع 


مکانی هدایت الکتریکی است. 

مطابق با نتایج مطالعه حاضر در مطالعه‌ای که توسط گومز و 
همکاران (2020 ,۵1 2 5ع00۳) برای پیش‌بینی محتوای 506 
انجام شد نتایج تایید کردند که مدل ۳ در هر پنج عمق مورد مطالعه 
نسبت به سایر الگوریتم‌های مورد بررسی شامل (کوبیست» 01/۷15 
بالاترین مقدار ضریب تعیین و کمترین خ 
مقایسه دو مدل "1 و 5۷5 برای پیش‌بینی 12 توسط وو و همکاران 
(2018 .۰۶1 ۷۷) نشان داد که مدل 1۳ با 182 بالاتر و 2۷۲51 کمتر 
دارای صحت بیشتری برای پیش‌بینی ۳ خاک است. برخلاف نتایج 
مطالعه حاضر که نشان می‌دهد پیش‌بینی ویژگی‌های ۳ و 80 
توسط هر چپار مدل یادگیری ماشین 16۳ 121 ۶۷۶ و ۸ در 
لایه زیر سطحی دارای صحت بیشتر و خطای کمتر است نتایج مطالعه 
دی آنتا و همکاران (2020 .۵1 6 ۸06۵ 16) نشان داد که پیش‌بینی 
0 توسط مدل ۳ در عمق ۲۰-۰ سانتی‌متری دارای صحت 


خطا را داشته است. 


بیشتری نسبت به عمق ۵۰-۰ سانتی‌متری بود. با همچنین نتایج 
مشابهی توسط نبی‌اللهی و همکاران (2017 ,۵ ۶ نطقاامن۲20) و 
ژانگ و همکاران (2017 ,.0 21 2202) گزارش شد. 
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جدول ۳- نتایج اعتبار سنجی پیش‌بینی مکانی شوری و کرین‌آلی خاک با اتفاده از مدل‌های یادگیری ماشین 
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بیشترین تغییرات 12 برای پیش‌بینی ۳ و 506 در هر دو لایه 
سطحی و زیر سطحی توسط مدل ۳ برابر با 7/۱ محاسبه شد که 
نشان‌دهنده تنوع کم در پیش بینی است. نتایج این مطالعه تائید می‌کند 
که مدل 16۳ یک الگوریتم یادگیری ماشین قدرتمند برای پیش‌بینی 
است که با نتایج ناوار و معز (2017 ,۱۷۲۵۵۵26 0طه ۱2/2۲) مطابقت 


دارد. 


تغییر پذیری مکانی 6 و 906 

حداقل و حداکثر ۳ در لایه سطحی به ترتیب ۰/۴۴ و ۸۲/۸۵ 
(شکل ۵-الف) و در لایه زیر سطحی به ترتیب ۰/۳۵ و ۷۷/۱۷ (شکل 
۵-ب) دسی زیمنس بر متر مشاهده گردید. در بخش‌های شمالی» شمال 
شرق و شمال غرب به سمت مرکز منطقه مورد مطالعه میزان ۲ کم 
بوده و بیشترین مقدار خطر و محدودیت شوری در مناطق جنوب و 
جنوب شرقی پیش‌بینی شد. بیشترین و کمترین مقدار 90067 پیش‌بینی 
شده در لایه سطحی ۰/۳ و ۲/۱ (شکل ۵-ج) و زیر سطحی به ترتیب 
پرابر با ۰/۱ و ۱/۳ (شکل ۵-د) درصد مشاهده گردید. با توجه به شکل 
۴-الف و ۵- الف و ۵-ب در منطقه مورد مطالعه بین تغییرات ارتفاعی 
و مقادیر ۳ پیش‌بینی شده روند عکسی وجود داشت. با توجه به تشریح 
واحدهای فیزیوگرافی در بخش مواد و روش‌ها قسمت شمال منطقه 
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ارتفاع بیشتر (دارای واحدهای فیزیوگرافی تپه و فلات) و در قسمت 
مرکز و جنوب ارتفاع کم است بنابراین در قسمت با ارتفاع کمتر میزان 
بیشتر است. به‌طوری که در ارتفاعات (۱۷۴۸ متر) مقدار ۲ تا 
مقادر کمتر از ۱ (15/00) و در ارتفاعات کم شامل اراضی پست با سطح 
آب زیرزمینی بالا (۱۱۳۹ متر) مقدار ۳ به بیش از ۸۰ (19/5) نیز 
می‌رسد که با مشاهدات میدانی و وضعیت واحدهای فیزیو گرافی منطقه 
منطبق می‌باشد. همچنین در قسمت شمال شرق» شرق و مرکز منطقه 
مورد مطالعه بیشترین محتوای ۹0۳ مشاهده شد که به علت وجود 
کاربری کشت آبیء دیم و همچنین وجود مراتع غیر شور می‌باشد در 
حالی که در قسمت جنوب و جنوب شرق, به علت وجود دشت‌های شور 
و عدم وجود پوشش گیاهی» 506 کاهش يافته است. به طور کلی 
0 در لایه سطحی تجمح‌یافته و با افزايش عمق خاک از میزان آن 
کاسته می‌شود که موسوی و همکاران (2022 ,.اه 61 ۷]۵۵5۵۷1) و فو 
و همکاران (2021 ,0 6 ۲۷) نیز نتایج مشابهی را گزارش کردند. به 
نظر می‌رسد در منطقه مورد مطالعه وجود رژیم‌های حرارتی ترمیک و 
رطوبتی آکوئیک سب برقرار شدن جریان موئینگی شده و املاح محلول 
از عمق به سطح به خصوص در فصول خشک سال انتقال می‌یابند 
(2021 ,.۵1 61 ۷۲۵۵5۵۷1 
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شکل ۵- توزیع مکانی الف: ۲ در لابه سطحی با استفاده از مدل 87 ب: »۳ در لایه زیر سطحی با استفاده از مدل 8۲ ج: »۹0 در لابه سطحی 
با استفاده از مدل 13۳ د: »50 در لایه زیر سطحی با استفاده از مدل 1۳۳ 
,۱0۵0061 1۳ 6ظ) عصزعت مرها آنمز طامای م6 رز ۲ :۵ وا0مصه ۱۳۲ مط عصنعت مرها آذمی رما م6 صنز ۳ جه ممتاتماتاعنل امتادمرک -5 مساعز۲ 
۵۵ ۱۱۳ م۲ امن ب۵رها ازمو ای ماه صز 90۵ 0 ملعم با عمط مهصنمته مرها آامز رما مص 90۵ 


نتیجه گیری 

این 7 تحقية با هدف بررسی تغییرات مکانی سطحم (۳۰-۰ سانتی - 
متری) و زیر سطحی (۶۰-۳۰ سانتی‌متری) ۳ و 506 توسط چهار 
مدل یادگیری ماشین (۲1۲ ,۸۲ ,5۷5 ,۳ در دشت قزوین 
انجام شد. از میان شش متغیر کمکی توپوگرافی و سنجش از دور مورد 
استفاده پارامترهای ۵۱۵7210108 ععبلنا و ۷۲۷۲۳۲ مهمترین 


متفیرهای پیش‌بینی کننده 6ظ و 506 در هر دو لایه سطحی و 
زیرسطحی بودند. که این نتیجه بیانگر این است که در منطقه مورد 
مطالیه هی کمک که فا کی ان کوییگرانی دز مشاه 
مطالعاتی هستند نسبت به شاخص‌های سنجش از دور از پیش‌ران‌های 
قوی‌تری برای پیش بینی این دو ویژگی مورد مطالعه (10 و 90) 
بودند. 

دو مدل یادگیری ماشین ۳ و :121 نسبت به مدل‌های ۸۳ و 


رضایی و همکاران. نقشه‌برداری رقومی تغییرات سطحی و زیرسطحی کربن آلی و شوری خاک در بخشی از اراضی دشت قزوین ... 


۴ تنایج دقیق‌تری را در تهیه نقشه پیش‌بینی مکانی ۳ و 5060 
مکانی ۳ نشان داد که بخش جنوب و جنوب شرق منطقه هم در 
سطح و هم در عمق دارای شوری بسیار زیاد بوده که محدودیت جدی 
برای کلیه کاربری‌های زراعی و باغی در منطقه مورد مطاله می‌باشد. 


۳۳۹ 


اراضی این منطقه برای کشت محصولات موجود در الگوی کشت 
منطقه نسبت به سایر بخش‌های منطقه باشد. در پایان از دیدگاه 
کاربردی نقشه‌های ۳ و ۹06 تهیه شده در اين مطالعه می‌توانند به 
عنوان اطلاعات پایه برای کاهش اثرات منفی شوری خاک و همچنین 
اتخاد تصمیمات مدیریتی برای افزايش بهره‌وری توان حاصلخیزی 
خاک‌ها در اراضی مورد مطالعه قرار گیرند. 
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ن بعصامم2 501 لماتعن 10۲ امما 2 قح 5۵۵00۲09600۲ ممصمامع۲۵1 1۲۳1161780 .(2009) .یه ررمصاه :۱۷۲03 ک ۵۰۷۰ ,م1۵55 
64) متا ۳/۱ عصتمممصه آزمی آمازع زر[ 

(۵11600۱-16607۷ ۷ کی و.گ ماک یگ بصوطاتتطمصمطه ریق ,۷۲۵052۷1 و.] .۷۲ متتقااهی بک بتهکصهنمصهه رگ رتصهطامطه 
0ص و10۳6 حممل‌صه ‏ فمماطمنمه اومتهمهص کر ممتقوم2ع۳ فقمعم ما0 مصتفتا مماه۲ممه7ن وم ۵0۵1182 :(2020) 
6 ۱( ۵ .66 ,(1)] 7 ,۸۸۱۱۵0۵/۱6 .28817515 ۵4۷6 مومع :وعصتطمقصظ ۷6۵۵۵۲ ٩۱۳۵۵0۵۲۲‏ 
۵ مولوها اتالامتنه7 المع ۵۶ )جمجصووع‌عفظ ,(2019) .ویک بقفقظ عک ۲۱۸ رتفهل پیایک مطتطقفطه یقت مطعع9۳11 
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رضایی و همکاران. نقشه‌برداری رقومی تغییرات سطحی و زیرسطحی کربن آلی و شوری خاک در بخشی از اراضی دشت قزوین ... 


905-۰ ,(9)01 ۲ ,۸۱۵۱۷۵5 ۳۱۵۲۱۶ ۵۱/۵۲ ر0115و 7121 10ج قمع بو منامتصطهعا از تاکزدع۳ 

11 ب,طاوه۲ز مه مموگ ۵۶ فعتنتهع۲ عملم۱۷۲20 2020(۰) ایو رمفتام] ی ربا بصقتطعهم۷ ,.۲ مصونلعفصنهو ریظ بتاقه ]1 
67-۰ و(۵)1 ل ۲۲۱۷۱۲۵۳۸۵۱۵ 6 ۱۱۱۵۱ ۰ ت۲۱ وه ۵۶ )و۲ طا ممباموتصطهع ۱ لمعم‌تاوتاهادمعن) مصتوتا ها مع۳۴۵111 
-560)1061 8۸۵ 09۵1-001271760 عصتویا مصتم مفصط واتصتاقی [رم9 ,(2019) بع بتتمطام]۲ ی ببا۷ بتاماصهوه۲۱ ,,۷۲۰1۷۲ رمع ]" 
.237-۰ «(40)1 ,967161118 1-01 "0 0۵ ۵۵۵ 11۳128261۳۷۰ 1 
7 0 143 1080/0 .105://001.0۲2/10ظ 

.(2021) ۲۰ ,001669 عک یه ,۱۷۲۵۹۵۷1 و.] مقصهتطظ ر.ظ رمصما! .لا رصقتصعجصهه0 1 ویک ر6تصصطم‌و ویک بللته[تط206-۷]۵ 12 مه ] 
٩0۵۵0۲ ۲‏ 24 ممتامصم]فصه م۷۵۷۵ عصلوهطا فحمتععم 20 طا واتصتاهی ام ۵۶ ممتامنلعم ۵202و معط عصتهمم‌جط] 
13 ۵5://001۰0۲82/10۰1016/1.826006۳۳98.2020) 114793 ,383 ,06006716۵ .38۵0615 ۲۵۵۲685109 

۳۲۲1۳۲1 ۵۶ همه مصتعه ممتاحصصتامی «اتصتامی ارم عموامواتای 0ع1آم‌صصزد ۸ 2021(۰) ۰۱۰ بتت/11۳ عک بیط بلطا۳10۵] 
رهاظ رتمعمصمرن؟ ما امه ومت عمط ۵۲ مممهاد باه ۶۵۲ بداعل ماتااماهه اهتاممموت‌انص ]9۳۲ 20 کف 
1/۵5۹://001.0۲2/10.0000.009 .3905-3919 ,(32)14 ,۵۷۵۱0۵۵ ۳۸۵8۵۲۵۵۵۱۵0۸۵ 1۲۵۳۱۵ :۲8012 

6۷۵۵۵01 ۱۱۵5 0۲۵۵ ۵06 ۸ ۲۷۵۸۵۵۵ 2006(۰) .۳۲ متا ی .۲.۷۷ ر108۵5 .نا بل۱۷2۱۵۳۵ ,.نا بطم۵112 ۷۷ 
0۴۱6۵1۱0۲۹۰ 0۵ 2۵۲۵۱۵۵200۲ 21105 

عصتطمقظ مصتوها مصتممفمط باتصتاهی ارم ,(2021) ۲۱ رمصعط2 ی ر.] ,۷۷202 .۷۷ بلالنا بب) مطالا وبا من ول ,۳۵92 و.ل ,۱۷۷22 
۱۹ ۱ اد ۸۱۸۱۱۹۱۱۱۱۱۱ ۱۰۵/۸۱۹۱۱۸۰۱ ۱۱۰۱9۱۱۱۱۱۰۹۰۱ 
05 و( 

نص .عرهب ای آخمی 0 ۱۱0۵600 اه م000۵ 000 ک ۷۵۲۱۵۵۵۵ 5۵/۵ 1985(۰) ۱:۳۰ ,عط۷۷۱۱1 
.166-4 .00 ,۲۷ ۷6۵۵۵5 فعض ,اتات2720 ۷ 502021 5011 

2016(۰) .6 بتالنآ ک و۷۰ مهعطه ,۲ .۵ متطوتهتام-ل۸ انط20ع۲ وب ۷۲۷۰ ,متکقط-۸1 ریک ۸ بلععصصتفطن ۱۷ بیی) بقمع0ا2 ,۷۷۰ ,۱۷۷۷ 
۵ ۱ بهلصطها0ومو ۱۷۲ آهصه ۲ صا ممتوومتوع: عصتصجع1 مصتطممصه ها مصتم‌مرمصظ فص ممتامتمنم ماتصتلمی تم 
10۰8 /005://101.0۲2 .4005-4014 و(29)11 ۵۷۵۱۵0۴۵۳۸۱۵۲ 6 ۵۵۲۵۵۵۱0۲۷( 

صح طا صمحاتقی متصقعه موز ۵ ممتام۳۳۵0۲ (2017) .) ,م6 ی لا رعطقط2 ویک رقصه۷ بیظ ما۷ بط بلا۷۷ بلط ,عطهط2 
حصمصه1۳ مص فص فصمنووم‌عع۱ بعمصناً مام‌ت انامه هم مموتنهم‌جومی ۸ تفص ممافمع ۵۴ عمصمم2 ومتمصهآ۲6۵ 0عمهصفظط ۷۵1۲ 1قصم 1۳ 
7۸6 010068۱۷ [/60۵5://01.۰0۲8۵/10.1016) 704-713۰ ,592 ۳۱۷۳۵۱۱۵۵ 101۵ 0۴1۵ 566۵۵ :20061 )8ع0۲] 

0 20201210100 20 ععقممموع۲ صقام ]۵ فحصوتصهطمع۱۷ ,(2020) .5 بقاهطاهطاه ع ریک بتاط2 بن) مقط0 با رعصقط2 بن) ,7120 
۰30188۰20207[ /605://01.۰0۲8/10.1016) 100017 :(1) ۲ ۲۱۱۳۵۵۲۷۵۸۱۵۷ 726 .ب)اتصتلهه 5011 

متصوعه آلمو ۵۴ موصتم‌ مق لم) تون جمتاتاآنعع۲۲۱2-۲6 2020(۰) ب۸۵۸ بطامعتاصضا ی رین ر0عقها۲ و.ل رطظ ول بصع .۷ بصع ,.] بتا20 
عصتصجع1 مصتطعمصه مم 0عوجه هم 2-امصتاجمو همع 1 -امصتاصمو رععتتاه۷ 0۵ ۱۳۱۷۲ مصنوها ممعمعنه لقاما انم صح ممماتنوم 
4 2020. 16/5610060۷ 10۰10 /6۵6://001۰0۲۵ :138244 72۵ 1۲۱۷۱۳۵۳۵۸۵۸۸ 10161 126 0۴ 567۱66 .21201110105 
متصحعه 9011 من 6۲6۵۲5 ععصقول م6مصصتآن 20 وبا 0صرض ,(2019) ۷۲۰ بتانا عک و2 ,28لا و2 بط بط یلصتججه) تقل ,۷.۰ :700 
1230-18۰ . ,64 ۱۳۳۵۱۵ 10101 ...126 0 56۵76۵ بصن .. افحعطات۱0 .. 0صه ... طاتملا.. صا.. صمهاوع 
۰16 7 1016/0110 10۰ /005://001۰012ظ 

عصتومهه مامصعن آهاممموتالنهه مصنوا مصتم‌مفصه ممطاعی متصععتن م5 ,(02019 ۲۰ ,220070۷ تک بط الط ,.ظ م7121 
۰ ,(24) ۲ ] ,5۵۳۵/8 ۲۵۱۵/۵ .006 ]0و6 تحامموو فصه لقتادمد من اه ماع ۶ اتاتطاد ممتام‌نلعظ تعاعل 
۰ و( 
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